Are black holes in alternative theories serious astrophysical candidates? 
The case for Einstein-Dilaton-Gauss-Bonnet black holes 
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It is generally accepted that Einstein's theory will get some as yet unknown corrections, possibly 
large in the strong field regime. An ideal place to look for these modifications is around the vicinities 
of compact objects such as black holes. Our case study here are Dilatonic Black Holes, which arise in 
the framework of Gauss-Bonnet couplings and one-loop corrected four-dimensional effective theory 
of heterotic superstrings at low energies. These are interesting objects as a prototype for alternative, 
yet well-behaved gravity theories: they evade the "no-hair" theorem of General Relativity but were 
proved to be stable against radial perturbations. 

We investigate the viability of these black holes as astrophysical objects and try to provide some 
means to distinguish them from black holes in General Relativity. We start by extending previous 
works and establishing the stability of these black holes against axial perturbations. We then look 
for solutions of the field equations describing slowly rotating black holes and study geodesic motion 
around this geometry. Depending on the values of mass, dilaton charge and angular momentum of 
the solution, one can have differences in the ISCO location and orbital frequency, relatively to black 
holes in General Relativity. In the most favorable cases the difference amount to a few percent. 
Given the current state-of-the-art, we discuss the difficulty to distinguish the correct theory of 
gravity from EM observations or even with gravitational wave detectors. 
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I. INTRODUCTION 



Einstein's theory of General Relativity (GR) has 
passed numerous consistency and experimental tests in 
a spectacular way jli]. Nevertheless, it is a general con- 
sensus that GR will get modified at some scale, if only 
because GR has resisted all attempts at its quantization. 
Moreover, other theories of gravity also pass the experi- 
mental tests, some with a better "quantum behavior" and 
should be taken as serious candidates. Unfortunately, 
the majority of these alternative theories are vastly more 
complex than GR in their full-fledged form. It is thus 
not surprising that progress in understanding the exact 
differences between one and the other and specially dif- 
ferences one can measure experimentally have been slow 
and mostly focusing on the weak-, far-field behavior. 

One of the candidates for a theory of quantum gravity 
is string theory [2]. Since it is still difficult to study ge- 
ometrical settings in superstring theories, most analyses 
have been performed by using low-energy string-inspired 
effective theories Q. Thus, quantum gravity predictions 
can be tested by studying modifications of GR due to 
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some low-energy truncation of string theory [J]. Typ- 
ically the effective theories are supergravities involving 
not only the metric but also a scalar field (the dilaton) 
as well as several gauge fields [5;] . One of such theories is 
the one-loop corrected four-dimensional effective theory 
of the heterotic superstrings at low energies and a sim- 
ple particular case is known as Einstein-Dilatonic-Gauss- 
Bonnet (EDGE) theory (see for instance Ref. and ref- 
erences cited therein for a nice and concise introduction 
to this theory). In EDGE theory the gauge fields are 
neglected and only the (spacetime-dependent) coupling 
between the dilaton and the gravity is considered, with 
the anomaly-cancelingChern-Simons term also neglected 
(see for instance 0,H,H,[l^ for work taking this term into 
account). At the first order in the Regge slope, a' , higher- 
derivative gravitational terms such as the Gauss-Bonnet 
(GB) curvature-squared term are present in the action, 
hence the name. The GB terms avoid some pathological 
features, for example the theory is ghost-free. Since the 
equations of motions are still of second order, EDGB the- 
ory provides one of the simplest consistent high-energy 
modifications to GR. Even though it does not seem to be 
a viable cosmological model [ll| , the Parametrized Post- 
Newtonian expansion of this theory is identical (to 
lowest order) to that of GR [H, , which means that it 
passes all Solar system-like experimental tests of gravity. 
Differences arise only from full nonlinear effects and the 
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ideal place to look for these is near compact objects such 
as black holes. 

Dilatonic Black Holes (DBHs) do exist in EDGE the- 
ory pisl . [l^ [iTj . They have a regular event horizon and 
the geometry is asymptotically flat at infinity (recently 
DBHs in higher dimensions [18] and with non-flat asymp- 
totical geometries 19] have also been studied). In what 
follows we refer to DBHs in EDGB theory only. DBHs 
are interesting for many theoretical and practical rea- 
sons. DBHs can evade the classical "no-scalar-hair" theo- 
rem , and be dressed with classical non-trivial dilaton 
hair. This is a direct consequence of the GB term [l^ 
and opens up the exciting possibility of ruling out such 
objects and theories by testing the no-hair theorem, ei- 
ther by gravitational-wave observations of ringdown f^T] 
or by observations of hi ghly eccentric orbits around su- 
permassive black holes [23]. The lack of experimental 
data and some problems on the theoretical side make it 
very difficult to explore string theory and other theories 
of quantum gravity. Since quantum gravity effects are ex- 
pected to play a significant role in strong gravity regime 
such as in cosmology [2^ and in black hole physics, the 
investigation of BHs in EDGB theory can shed new light 
on some aspects of quantum gravity and/or be used to 
develop testable predictions of the theory. 



A. Astrophysical implications 

Our purpose here is to begin exploring differences be- 
tween GR and alternative theories that could be exper- 
imentally tested through astrophysical observations, in 
the strong field regime. Thus, even though the theory 
with which we will work is only a first order trunca- 
tion of the full action, we will elevate it to the status 
of a complete theory, therefore we place no restriction 
on the Regge slope. The first important problem con- 
cerns the stability of DBHs. Credible alternatives to 
the Schwarzschild and Kerr metric of GR must be stable 
spacetimes. It has been shown [U,!!^ that DBHs are sta- 
ble against a small subset of all possible perturbations, 
linear radial perturbations. In this work we character- 
ize completely half of the possible degrees of freedom, 
by studying general axial perturbations. We find that 
DBHs are stable also against these perturbations. By 
itself this is an interesting result confirming that high- 
energy contributions lead to viable alternatives to classi- 
cal BHs arising in GR. The viability of DBHs poses the 
following question: can one devise observational tests to 
discern a DBH from a classical BH? In classical Einstein- 
Maxwell theory BHs are characterized by three parame- 
ters [l^: mass M, electric charge Q and angular momen- 
tum J = aM ^ M^. Astrophysical BHs are likely to be 
electrically neutral because of the effect of surrounding 
plasma [27[ , and therefore tests of alternative theories of 
gravity can in general focus on rotation alone. The task 
is still highly non-trivial and as we mentioned earlier, 
strong-field effects must be searched for in theories that 



are not already ruled out by Solar-system experimental 
data di]. 

Most if not all of present-day astrophysical observa- 
tions related to compact objects, concern directly or in- 
directly the motion of matter. Thus, a study of geodesic 
motion around compact objects in alternative theories is 
of utmost importance. Geodesies convey very important 
information on the background geometry. In particular 
circular orbits whose radius is close to the horizon may be 
extremely useful, because they already probe strong field 
regions. They can be used to compute the "innermost- 
stable-circular-orbit" (ISCO), a notion which is very im- 
portant for interpretation of the experimental data con- 
cerning astrophysical black holes. For the Schwarzschild 
spacetime, nsco = 6M, while for an extremal rotating 
Kerr geometry, risco = 9M for co- and counter- 
rotating circular orbits respectively. Measurements of 
the ISCO are also useful to evaluate the angular mo- 
mentum of Kerr BHs (28j . Current methods to measure 
the ISCO include spectral fitting, quasi-periodic oscilla- 
tions and relativistic iron line measurements [28i]. Here 
we show that differences in the ISCO of slowly rotating 
DBHs and classical Kerr black holes can be significant, 
depending on the coupling parameter. These differences, 
which may be detectable in near-future experiments, are 
likely to increase for highly spinning black holes (which 
are unfortunately out of the scope of this work). In fact, 
it seems an exciting possibility that the GRAVITY ex- 
periment [2^ . designed to make precision measurements 
of orbits of stars in the neighborhood of the black-hole 
in the center of our galaxy, might already be able to 
discriminate between these black hole solutions and GR 
solutions, or otherwise impose stringent bounds on the 
coupling parameter. 

Current techniques to evaluate the angular momentum 
of astrophysical compact objects are based only on ISCO 
measurements 28]. One needs to assume a particular 
theory of gravity in order to evaluate J. Thus discerning 
a DBH from a Kerr black hole by ISCO measurements 
is not an easy task. Future gravitational wave experi- 
ments will provide a viable method to measure M and J 
independently [s^]. The analysis below suggests that a 
possible deviation from the expected ISCO in GR can be 
explained in term of dilatonic charged BH. Therefore in a 
near future, gravitational wave astronomy may offer a the 
possiblity to explore string theory-inspired modifications 
of GR. 

Finally, null unstable geodesies are closely related to 
the appearance of compact objects to external observers 
[31] . [3^ and have been associated with the character- 
istic, or quasinormal modes (QNMs) of BHs [H, [s^ ]. 
Quasinormal modes are very important in devising ex- 
perimental tests of GR and for gravitational wave as- 
tronomy. Measuring QNM frequencies may provide a 
definitive proof of the existence of BHs in GR and it 
could be useful to study corrections to GR too. Thus the 
analysis of geodesic motion around a DBH can shed new 
light on various and important aspects of high-curvature 
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corrections to gravity. 

Other theories might also suggest a different ISCO lo- 
cation, different quasinormal modes, etc. Thus, a devia- 
tion in these quantities is not a verification of a particular 
theory, but is a first step in understanding what different 
theories and scenarios predict in the strong field regime, 
which could potentially discriminate GR from other al- 
ternatives. 

The paper is organized as follows. In Section |TT] we 
briefly review the main aspects of BHs in EDGB the- 
ory. In Section IIIII we prove the linear stability of DBHs 
against axial perturbations. We proceed by studying 
slowly rotating DBHs in Section HVl We prove that such 
slowly rotating solutions do exist and we characterize 
them, including a discussion on the ergoregion in these 
spacetinies. Section |V] discusses geodesies in both spheri- 
cally symmetric and slowly rotating DBHs as well as the 
possible experimental tests which can be performed to 
discern a DBH from a Kerr BH. We compute the ISCO 
dependence on the angular momentum and the QNM 
frequencies for a spherical symmetric DBH in the eikonal 
limit. Conclusions are discussed in Section lVTl In our no- 
tation, we use the signature (H ) and the curvature 

tensor is defined by R^j^. 



Here, 



II. SPHERICALLY SYMMETRIC BHS IN 
EINSTEIN-DILATON-GAUSS-BONNET THEORY 

We consider the following low-energy effective action 
for the heterotic string [3] 



5 = ; / d^x. 
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where 



n 
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(2.1) 



(2.2) 



is the GB invariant, a' is the Regge slope and g"^ is some 
gauge coupling constant. We set 5 = 1 for the rest 
of this paper. String-inspired 0{a') corrections to Ein- 
stein's gravitation are included in the action (12. while 
gauge fields and matter are omitted for simplicity. We 
also note that there is some arbitrariness in the coup ling , 
depending on which frame we take as fundamental [l8|, 
we choose to keep the e*^ coupling in line with previous 
works. We shall refer to the action (|2.ip as EDGB theory. 

The dilaton field and Einstein's equations derived from 
(1^ are 
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GB ' 



(2.3) 



Gp, = ^d^cl)d,<l> ~ ^g^.(9p0)' - a'lC^, , (2.4) 
where G^^ — R^^ — ^g^^R is the Einstein tensor and 

= [g^^pgux + 9^.x9.pW^^''^yg[R''^^0^J] . (2.5) 



i-g)- 



Oijk _ _ 



^ kX — 'I ^panX , ~ 8 ' 

From the right-hand-side of the modified Einstein's equa- 
tion (|2.4[) , one can construct a conserved "energy momen- 
tum tensor", V„r^'' = 0, 



(2.6) 



In Ref.[15] is shown that the time component of —T^^^ 
which in Einstein's gravity would correspond to the local 
energy density £, can he negative. The reason is that, 
as a result of the higher derivative GB terms, there are 
contributions of the gravitational field itself to T^^. The 
positiveness of — Tqo is one of the main assumptions of 
the no-scalar-hair theorem poj which can be (and indeed 
is) evaded in EDGB theory. 

We now focus on BH solutions in EDGB theory, con- 
sidering the following spherically symmetric ansatz for 
the metric 



r^{d0^ +sm^ 9 dip^). (2.7) 



The equations of motion derived from l|2.3p and from 
(|2.4p can be found in Ref. In a static, asymptotically 
flat geometry, black hole solutions exist only if [l^ 
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where and (f)h are quantities evaluated on the hori- 
zon. In particular black hole solutions may exist only for 
a' > 0. From the asymptotic behavior of the fields one 
can extract the ADM mass, M and the charge D. As 
shown in Ref. M and D are not independent quanti- 
ties, thereby leading to the secondary nature of the dila- 
ton hair [3a|. These black hole solutions are uniquely 
characterized by two parameters {(j)h,rh), which corre- 
spond to a unique choice of {M,D). The equations of 
motion remain invariant under a shift 4) c/) + (po and 



a simultaneous radial rescaling r 



^ji>oli_ ^ conse- 



quence of the radial rescaling, the two other asymptotic 
parameters, M and D, are also rescaled according to the 
rule M and D De'^"/^. Due to the above 

invariance it is sufficient to vary only one of rh and (f>h- 
Following we choose to keep rh fixed and to vary (j)h- 
Typical background fields are shown in Fig. [T] Differ- 
ences in the metric coefficients occur only very close to 
the horizon. We checked our numerical solutions repro- 
ducing results shown in the Table 1 of Ref. [lB| ■ 

After the rescaling, equation (|2.8p can be written in 
terms of the coupling constant 



< < 0.691. 



(2.9) 



4 



The maximum value a' /M"^ ^ 0.691 corresponds to 
D/M ~ 0.572. For larger values of the coupling con- 
stant no DBH solution exists. The dilaton charge, Z?, as 
a function of q'/M^ is shown in Figs. [21 where it is com- 
pared with the analytical solution in the a' — > limit, 
D/M = a'/{2A'P) [ai. 




g„(r) EDGB 
) EDGB 
&.(r)GR 
^;,T(r) GR 



1.4 1.6 
r/rh 



FIG. 1: Metric components gtt and grr for a DBH (solid lines) 
compared with an equal mass Schwarzschild hole (dotted 
lines) (a'/M^ ~ 0.691, which corresponds to D/M ~ 0.572). 




FIG. 2: The dilaton charge, D, as a function of (re-scaled) 
a'/M^ (solid line). The whole range of a'/M^ in which a 
black hole solution exists is shown. Numerical solution be- 
haves as D/M ~ a'/{2M^) in the a' -* limit ;3e] (dashed 
line) . 



(with no dilaton), and generalized by Moura and Schi- 
appa 61 in the context of Riemann tensor, TZ^, correc- 
tions. Perturbations of four-dimensional DBHs were con- 
sidered by Kanti et al [2?] and Torii and K. i. Maeda [25j . 
Unfortunately these authors considered only a very spe- 
cific type of perturbations, here we want to generalize 
their results. Consider therefore perturbing the space- 
time in the following linear way 



AO) 



(3.1) 



where e ^ 1, gl^J and ^0(2^'') are the background fields, 
while and S(j){x'') are the perturbations. The back- 
ground metric gjj^J and dilaton field (I)q{x'') are given by 
the numerical static black hole solution described above 

M. 



A. General formalism 

To study this problem we use the approach first de- 
scribed by Regge and Wheeler [i^. After a decompo- 
sition in tensorial spherical harmonics (4ll . |43 | , the per- 
turbations fall into two distinct classes: axial (odd) with 
parity (—1)'^^ and polar (even) with parity (—1)', where 
I is the angular momentum of the particular mode. The 
theory described by (|2.ip is invariant under diffcomor- 
phisms"'^, as Einstein's theory is. We can then use the 
gauge freedom in order to simplify the elements h^^^. In 
the classical Regge- Wheeler gauge the canonical form for 
the metric perturbations is (see also Ref. ^43j): 
- axial perturbations: 



/lo 

/ii 



ho hi 



(sindde) Pi, (3.2) 



polar perturbations: 



Hi 








r^K 








r^A'sin^ 6 



'Pi. 



Where Pi = Pi{cos9) is the Legendre polynomial with 
angular momentum I and ho, hi. Hi, H2 and K arc un- 
known radial functions. 

Perturbations of the dilaton field, S(f), do not appear 
in the axial equations (see also The linear stability 

analysis proceeds by mapping the system of the equations 



III. LINEAR STABILITY ANALYSIS 



The framework for a linear stability analysis of black 
holes in theories with Gauss-Bonnet terms was laid down 
by Dotti and Gleiser [s^, [H, [s^ in higher dimensions 



^ Despite the covariant derivative in eq. 1)2. 5| l , equations of motion 
do not contain higher derivatives of the metric because of 
the GB term (see equation (4) in [5^). This allows one to use 
the same gauge transformations first proposed in Ref. |40|| . 
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of motion for the perturbation fields under consideration 
to a stationary one-dimensional Schrodinger problem, in 
an appropriate potential well, in which the 'squared fre- 
quencies' uP' are the "energy eigenvalues". Instabilities, 
then, correspond to bound states, i.e. to negative en- 
ergy eigenstates or equivalently to frequencies lo with a 
positive imaginary component. 

Presumably because polar perturbations are extremely 
complex to analyse, Kanti et al. [23| focused on a certain 
subset, the radial perturbations. They found that the 
spacetime was stable against radi al p erturbations, their 
results being confirmed in Ref. [2^. Thus, there are 
good indications that the spacetime is stable under polar 
perturbations. We thus focus here on the other set of 
perturbations, axial perturbations, which as far as we 
know are not dealt with in the literature. 



B. Axial perturbations 

For axial perturbations only 3 nontrivial Einstein equa- 
tions can be obtained by substituting ()3.2|) and ()3.1|1 into 
(|2.4p . The zeroth-order equations are identically zero, 
due to the background solution, while the equations for 
the perturbations, /lo and /ii, read: 



hoir) + Al hi{r) + A2 h[ir) = , (3.3) 

h'oir)^^hoir) + l,lhi{r)^0, (3.4) 

Ci,iho{r) + C2jhi{r) + C3,ih'a{r) + 
Ci,ih[{r) + C5jK{r)^0, (3.5) 



where Ai, Bi^i and Cu depend on the radial back- 
ground function F, F', F", A, A', (/)o, 0o and </)g found in 
[iSj . Their explicit form can be found in Appendix [XI 
We observed numerically that equation (|3.5p is automat- 
ically satisfied as a consequence of the other two equa- 
tions and of background solutions. So we are left with 
a system of two ODEs for two unknown functions hQ{r) 
and hi{r). Eliminating from the first order equation 
p.3p we obtain a second order differential equation for 
hi which can be recast (see Appendix for details) in 
the following Schrodinger-like equation. 



u"{r) + [V{r) + uj^K{r)\ u{r) = , 

These functions are shown in Fig. [31 

The asymptotic behavior of equation 
pendix YK\ for details) is 



u" (r) 



(r - ThY 
u"{r) +Lo'^u{r) =0, 



u{r) = 0, 



(3.6) 

(see Ap- 

Th , (3.7) 
00. (3.8) 



The asymptotic solutions are 

u{r) - Mo(r-rh)3^Vi-^'»-'^'-^'> , ^ ^ , (3.9) 
u{r) ~ Uooe^"^\ r^oo. (3.10) 



Figure [4| shows the coefficients Vh and Kh as functions 
of a'/M^. For a Schwarzschild BH, T4 = 1/4 and K,, = 
4:AP. In the EDGB case one finds one finds Vh ~ 1/4 and 
Kh < in the whole range (j2.9p . Thus the asymptotic 
solution (|3.9p simplifies 



u{r) ^ UQ{r ~ Th) 



Th 



(3.11) 



Since lo is complex, the sign in equations (|3.10p and 
p.ll[) above has to be chosen so that the solution is reg- 
ular on the horizon and at infinity. For unstable modes 
Im(a;) > 0, thus the choice of the minus and the plus 
sign in eqs. (|3.10p and (|3.1ip respectively guarantees 
that the corresponding eigenfunctions will vanish at in- 
finity and at the horizon. We have searched for unstable 
modes using these boundary conditions. We integrated 
equation (|3.6p outward starting from the horizon until we 
found the eigenfrequency corresponding to a vanish field 
at infinity. We used a Runge-Kutta 4th order method, 
considering rh — 1 and different 4>h values corresponding 
to different BH solutions |15|. We systematically span 
the whole range (|2.9p for / = 2, 3, 4. We also randomly 
span other values for a' and I to no avail: no unstable 
modes were found. We checked the numerical accuracy 
of our results by changing numerical parameters such as 
Too. Our results strongly suggest that EDGB black holes 
are stable against axial perturbations. This completes 
the previous linear stability analysis 0, [2^ performed 
for radial perturbations. 



It would be of great interest to understand also polar 
perturbations and reduce them to a single master ordi- 
nary differential equation. Due to the complexity of the 
equations involved we did not perform such an analy- 
sis. We note that the present analysis is only possible 
because axial and polar perturbations decouple in tjiis 
theory, unlike for instance in Chcrn-Simons theory [4^ . 



IV. SLOWLY ROTATING BHS IN EDGB 
THEORY 

Astrophysical black holes are expected to be highly 
spinning because of accretion and merger events [H, |43, 
[Isl . l49j . Unfortunately the GB term is very hard to deal 
with for generic axisymmetric metrics and thus general 
rotating black holes are difficult to study without a full 
numerical integration of Einstein's equations. Expansion 
in a' is another promising approach, which has been suc- 
cessfully implemented for non- rotating DBHs 36]. Here 
we use another method, searching for solutions describing 
rotating BHs at every order in a', but in a slow rotation 
approximation. This method allow us to prove that such 
(slowly) rotating BHs do exist in EDGB theory. This 
strongly suggests that, in general, rotating DBHs exist. 
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FIG. 3: Top Panel: The potential V{r) for the axial per- 
turbations for a EDGB black hole. Different values of the 
angular momentum / are shown for q'/A/^ ~ 0.691, which 
corresponds to D/M ~ 0.572. Bottom Panel: The function 
K{r) for the axial perturbations equation for a EDGB black 
hole for a'/M^ ~ 0.691. 



In alternative high energy modifications of GR, such as 
Chern-Simons modified gravity BHs, rotating black hole 
solutions are harder to find [50, Isil. 

We follow Hartle's approach {52] which is based on 
axisymmetric perturbations of a spherically symmetric 
equilibrium solution and on an expansion of the per- 
turbed quantities in a power series of the angular velocity 
51 for frame dragging. The equilibrium solution is the one 
developed in Ref. [15| and explored in the previous sec- 
tions. To first order in Q the perturbed metric is 

(4.1) 

where = ft{r, 9) is the angular velocity d(j)/dt of an 
observer at (r, 9) freely falling from infinity. Since we 
consider only axisymmetric perturbations, we neglected 
the (j) dependence and we can expand 51 (r, 9) in terms of 
Legendre polynomials 



n{r,9) 



1=1 



1 dPi{cos9) 



smf 



d9 



(4.2) 



Both the metric and dilaton perturbations are functions 
of even powers of the angular velocity fl, thus we can 
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FIG. 4: The coefficients Vh (solid line) and Kh (dashed line) 
as functions of a' /M"^ in the whole range (|2.9|) . As for a 
Schwarzschild BH, Vh = 1/4, whereas Kh is always smaller 
than its Schwarzschild counterpart, Kh < 4M^. 



use the unperturbed metric and the unperturbed dila- 
ton field to find il-order corrections neglecting fi^-order 
terms. Using the expansion ()4.2|) only the {t, (^} compo- 
nent of modified Einstein equations (|2.4p is first order in 
O and, interestingly enough, it gives a separable differ- 
ential equation as follows 



where 



G2(r) 



n', 



G,{r){2-l{l + l)) 



0, 

(4.3) 



Gi(r) 
G2{r) 



2e^ + a'e* {A'(f>' - 2(0'^ + cj,")) , 

-eV(-8 + r(A' + r')) + 

a'e"^ (0'(6 - r(3A' - 20' + F')) + 2r</.") 

2r^e^ - 2a're'^</)' , 



(4.4) 



and A, (f), T are zeroth-order in fl, therefore given by 
the spherically symmetric background [isj . For a' = 
equation (|4.3p is equivalent to the standard GR form 
(equation (43) in Ref. [s^l with vanishing stress tensor). 

To solve the above equation we must specify the 
boundary condition for large r. Since the scalar field van- 
ishes asymptotically, the solution should approach flat 
space solution. Thus for large r we can define the angu- 
lar momentum of the BH via the following behavior, 



2J_ 



(4.5) 



It is worth to note that the dilaton field can introduce 
order corrections to the off-diagonal metric coeffi- 
cients. In the slowly rotating approach we use, these 
corrections are neglected. The asymptotic behavior of 
equation (14. 3p is 



2-1(1 + 1) 



(4.6) 
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whose general solution is 

ni(r)^ar-^-'+l3r-^+\ (4.7) 

where a and (3 are constant. With the asymptotic behav- 
ior (|4.5p . the equation above implies that only the I — 1 
mode is allowed. Thus the equation to solve is 

+ ^1^1=0. (4.8) 

The solution of the above equation is 

n^{r) = Ci +C2 Tdte^^'-.'^^'^^^^^^/'^-^^^" , (4.9) 

where the constants Ci and C2 are fixed asking for the 
asymptotic condition at infinity, eq. (|4.5|) . The analo- 
gous procedure to find slowly rotating Kerr solutions in 
classical GR is presented in Appendix |B] for complete- 
ness, which is basically a reproduction of the results of 
Hartle [s^l- There we show that this procedure leads to 
the Kerr metric in lowest order and furthermore that the 
I = 1 term is the only possible term in the expansion. 
Applying the same procedure to a static boson star [s^ . 
one c an p rove that such slowly rotating solutions do not 
exist dl]. 

In what follows we shall drop the index, = ili. The 
angular velocity ^{r) for different slowly rotating BHs is 
shown in Figs. ([5]). In the limit a' 0, one recovers 
the results for a small rotating BH in General Relativity, 
obtained perturbing a Schwarzschild BH, as described in 
the section below. Interestingly, the angular velocity of 
a BH in the EDGB can be ^ 40% larger than the one 
for a slowly rotating Kerr BH (see the bottom panel in 
Fig. ([5])) with the same angular momentum. A numerical 
relation which can be computed from Fig.(l5|) is 

M^n= Mn{rh)^Q.il^, (4.10) 

while for a slowly rotating Kerr BH the above relation 
is MUh = 0.25 J/M2. Thus, a DBH is more "compact" 
than a BH with same mass and angular momentum. 




FIG. 5: Top Panel: Angular velocity f2 as a function of r for 
different slowly rotating BHs in EDGB theory. From top to 
bottom: decreasing values for ol jM^ . For a' ^ or D — > 
solutions approach the slowly rotating Kerr BH (solid line). 
Bottom Panel: Angular velocity of the BH, Q.h ~ fl{rh), for 
a DBH (dashed line) with a' /M^ ~ 0.691 compared to the 
same for a slowly rotating Kerr BH (solid line) and a Kerr 
BH (dotted line). Up to J/M^ ~ 0.5 results for the slowly 
rotating Kerr BH reproduce the exact ones. The difference in 
the angular velocity between slowly rotating DBHs and slowly 
rotating Kerr BHs and can be as large as ~ 40%. 

Eq. (|4.1ip is topologically a torus. In the equatorial plane 
we have 

rn{r) = Ve^(^) . (4.12) 



A. Ergoregion and superradiance 

Ergoregions can develop in rotating spacetimes. The 
ergoregion is found by computing the surface on which 
gtt vanishes. An approximate equation to the location of 
the ergoregion [55] is 

= -e^'-'''^ +Q'^{rysm'^9, (4.11) 

which is expected to be a good approximation specially 
for very compact objects, such as BHs. The solution of 



The existence and the boundaries of the ergoregions can 
be computed from the above equations. The ergoregion 
width for a DBH compared to both a slowly rotating 
and a full rotating Kerr BH is shown in Figure [H] for 
different J values. The ergoregion width, W, for a Kerr 
BH on the equatorial plane is W/M = 1 - a/1 - {a/My. 
The ergoregion width for a DBH can be ~ 50% larger 
than the width for a slowly rotating Kerr BH. Therefore 
BH superradiance in EDGB is expected to be stronger 
than in GR. As expected the difference between slowly 
rotating DBHs and slowly rotating Kerr BHs tends to 
zero in the a' — > limit. 
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FIG. 6: Ergoregion width, for a DBH (solid line), a Kerr 
BH (dashed line) and a slowly rotating Kerr BH (dotted line) 
as function of the angular momentum J for q'/M^ ~ 0.691. 
Up to J/M^ ~ 0.3 results for the slowly rotating Kerr BH 
reproduce the exact ones. The ergoregion width for a rotating 
DBH can be ~ 50% larger than the one for a Kerr BH. 



V. BH GEODESICS IN EDGE THEORY 

In this section geodesies in the exterior of both non- 
rotating and slowly rotating EDGB black holes are dis- 
cussed. If we neglect back-reaction effects, which we do 
in the following, the geodesies in this spacetime corre- 
spond to paths followed by time-like or null particles. In 
principle back-reaction effects should be important for 
large bodies or strong dilaton fields. However, for small 
point particles around DBH black holes they should be 
negligible. From the analysis of the geodetic motion we 
compute orbital frequencies related to the ISCO for both 
non-rotating and slowly rotating DBHs and the QNM 
frequencies for spherically symmetric BHs in EDGB the- 
ory in the eikonal limit. These quantities can be directly 
measured and can be used as a promising tool to study 
high-energy modifications of GR. 

The motion of test particles around a DBH is deeply 
related to the coupling between matter and the dilaton. 
However it turns out that the coupling depends on the 
choice between two different frames, namely the string 
frame and the Einstein frame (for a discussion in scalar- 
tensor theories see Ref. [gII and references therein). The 
two frames are inequivalent. In the string frame the cou- 
pling between matter and the dilaton field is minimal, 
but Newton constant, Gn, depends on the coordinates 
and light bending is not correctly reproduced [6^ . Con- 
versely, in the Einstein frame the light bending is cor- 
rectly reproduced, since Gat is constant, but particles 
non-minimally couple to the dilaton field. This leads to 
a violation of the equivalence principle. This violation 
is small and compatible with the available tests of the 
equivalence principle [68| and it is indeed regarded as an 
important test in discerning the two frames. Which of 
the two frames is the physically relevant is an open ques- 
tion, which will eventually be settled by experiments. 



Following Refs. and [l6j we always assume the 

Einstein frame. In this frame the total action reads 



S'tot 



S -\- Sc + Sts 



(5.1) 



= S + Sc-m J dte'^^^-g^.- 



'dt~dt' 



where S is the EDGB action, eq. (|2.ip . m is the mass 
of the test particle and Sc is the action which describes 
some coupling between the test particle and the dilaton 
field. The constant b in the equation above is the coupling 
between the matter and the dilaton field. Its particular 
value depends on the specific theory from which the low 
energy theory comes. For low energy modifications from 
heterotic string theory, b = 1/2. We shall discuss the 
motion of a test particle keeping a general value of b 
and we only specialize to 6 = 1/2 when we discussing 
numerical results. As already mentioned, in the equation 
above we neglect gauge fields, such as Maxwell fields. 
Furthermore we neglect any coupling between the dilaton 
field and the test particle, S'c=0, as well as any back- 
reaction of the background fields. 

The string frame and the Einstein frame are related by 
a Weyl transformation (t^]. Since null geodesies equa- 
tions are Weyl invariant, the motion of massless particles 
is the same in both frames and it is described by the 
standard geodesies equations. 



A. Geodesies: non-rotating case 

Following Chandrasekhar [SB^ and Ref. [33| we consider 
a four-dimensional stationary, spherically symmetric line 
element 



ds^ = f{r)dt^ 



1 



dr^-r'^{d9^ + 



sm 



'■)■ (5.2) 



We also take a spherically symmetric dilaton and set 
H{r) = 260(r). Due to the symmetry of the above back- 
ground fields, the trajectory of a particle is planar, say 
9 = t:/2 and the following conserved quantities can be 
defined: the (dimensionless) specific energy of the test 
particle, E — e/m and the specific angular momentum, 
L = £/{m), where e, £ and m are the energy, the angular 
momentum and the mass of the test particle respectively. 
Let us restrict attention to circular orbits, for which the 
Lagrangian is 



2C = e^^'-) ( f{r) P - -^r' - r^^A 



(5.3) 



From the Lagrangian above, the equations of motion for 
the coordinates — (t, r, 6, (f>) read 



V(r) 



9ir) 



a2H(r) 



fir) 



E 
W) 



7^2 

X H(r) 
-if(r) 



, (5.4) 

(5.5) 
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where Si = 0, 1 for light-like and timelike geodesies re- 
spectively and the derivative is intended to respect with 
the proper time. 



1. Time-like geodesies 

We set Si = 1 in equation (15. 4p . The requirement for 
a circular orbit at r = rt is V{rt) — V'{rt) — 0, thus 



We choose the notation ft = f{rt), with the subscript t 
standing for timelike. Since the energy must be real we 
require 



2 + nHi 



> 0. 



2/t - nfl 

The condition for a stable circular orbit is 



(5.6) 



V' 



Vt - Ttfl 



f t r „ll 



{H[f] < 



(5.7) 




From equations (|5.5p we define the orbital angular veloc- 
ity 



(5.8) 



Equations above reduce to the usual ones for H{r) = 0. 

For the well known Schwarzschild case we have /(r) = 
g{r) ~ \ — 2M/r and H{r) = 0. In this case the condition 
(|5.6[) reads rt > 3M and from condition (|5.7p for stable 
orbits, we find rt > QM. For 3M < rt < 6M only un- 
stable orbits exist, hence the radius r = risco = 6Af 
is known as "innermost-stable-circular-orbit" (ISCO). 
The orbital angular velocity at the ISCO is Mfiisco = 
1/(6V6). 

To compute the above quantities in EDGB theory we 



set fir) = e^W, g{r) 



and H{r) 



where 



r(r), A(r) and (f>(r) represent the spherically symmetric 
BH solution found in Ref. [l^. We also specialize to the 
case b — 1/2. Numerical results for a EDGB non-rotating 
BH are shown in Table [T] and Fig. [T] The comparison 
with Schwarzschild BH (with the same mass) and inclu- 
sion of null geodesies is discussed in the section below. 
For completeness we show results both for 6=1/2 and 
6 = 0. The effect of a non- vanishing coupling is leading 
and it makes the angular frequency for a DBH smaller 



than in the Schwarzschild case (conversely the ISCO is 
larger). We find that the difference between the orbital 
frequency in GR and in EDGB theory can be as large 
as ~ 20%, depending on the coupling constant a'/Af^. 
In the a' ^ limit Schwarzschild results are recovered. 
The largest deviations occur when relation ()2.8p is sat- 
urated by chosing the appropriate (j)h- This corresponds 
to the maximum value a' /M'^ ~ 0.691, or equivalently 
to Z3/A/ ~ 0.572. It is worth noticing that the same 
occurs for a Reissner-Nordstrom BH, for which we have 
a relation between the electrical charge and the mass, 
Q/M < 1. Qualitatively, taking in account the coupling 
6, in EDGB theory the ISCO is always larger than it 
is for a Schwarzschild BH and the orbital frequency for 
a timelike geodesic is always smaller. Results for null 
geodesies do not depend on the coupling 6, as explained 
in the section below. In Table|T]we also show the quantity 



V = 



-E^isco 



So 



which is the binding energy per unit rest-mass at the 
ISCO. Because the accretion inside the ISCO is assumed 
to be in free-fall, this is equal to the integrated luminosity 
per unit rest-mass accreted, i.e. the radiative efficiency 
of accretion. This efficiency ranges from ^ 5.7% for a 
Schwarzschild black hole hole to ~ 42% for an extremal 
Kerr hole. 

Finally, we searched for static equilibrium solutions 
(following Maki and Shiraishi [7l[), i.e., a point particle 
at rest at a distance say, r = rg. For this, we set L = 
in Eq. EH and ask for V{ro) = V'{ro) = 0. We find no 
solution to these equations, meaning such a configuration 
does not seem to be possible for this theory. 



1.00 




FIG. 7: Angular frequency for circular timelike, fiisco, for 
6 = (solid line) and for b — 1/2 (dashed line) and 
for null-like, flc (dotted line) geodesies, normalized to the 
Schwarzschild value. 
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TABLE I: Some geodesies-related quantities for a spherically symmetric EDGB BH: orbital frequencies for timelike (fiisco), 
and null (f2c) circular orbits and the radiative efficiency, 77. We compare to the Schwarzschild case, for different values of a' /M"^ 
covering whole spectrum < a' /M'^ < 0.691. We also give the percentage (%) difference between EDGB and GR and we show 
results both for the coupling 6=1/2 and & = 0. Null-like geodesies do not depend on the coupling h. 





3 X 10"' 
3 X lO"' 
0.027 
0.238 
0.417 
0.545 
0.630 
0.677 
0.691 



D/M 





10~^ 
2 X 10~ 
0.013 
0.129 
0.246 
0.350 
0.441 
0.516 
0.572 



MQ.1SC0 [h = 0) (±%) 



0.0681 
0.0681 
0.0681 
0.0681 
0.0681 
0.0684 
0.0687 
0.0691 
0.0695 
0.0697 



6V6 

(~ 0%) 
(~ 0%) 
(~ 0%) 
(-FO.2%) 
(-hO.5%) 
(+1.0%) 
(+1.6%) 
(+2.1%) 
(+2.5%) 



MQ. 



ISCO 



(±%) 



0.0681 - 
0.0681 (■ 
0.0680 (- 
0.0677 (- 
0.0646 (- 
0.0617 (- 
0.0591 (- 
0.0569 (- 
0.0552 (- 
0.0539 (- 



6^ 

-0%) 

-0.1%) 

-0.5%) 

-5.0%) 

-9.4%) 

13.1%) 

16.3%) 

18.9%) 

20.9%) 



Kmc (±%) 



0.1925 ~ ^ 
0.1925 (~ 0%) 
0.1925 (~ 0%) 
0.1925 (~ 0%) 
0.1926 (+0.1%) 
0.1930 (+0.3%) 
0.1936 (+0.6%) 
0.1942 (+0.9%) 
0.1948 (+1.2%) 
0.1953 (+1.5%) 



7; (b = 0) (±%) 



3-2^2 



0.0572 - 
0.0572 (~ 0%) 
0.0572 (~ 0%) 
0.0572 (~ 0%) 
0.0572 (+0.1%) 
0.0574 (+0.3%) 
0.0575 (+0.6%) 
0.0577 (+0.9%) 
0.0579 (+1.2%) 
0.0580 (+1.3%) 



n (±%) 



0.0572 ~ 
0.0572 (- 
0.0571 (- 
0.0567 (- 
0.0525 (- 
0.0484 (- 
0.0449 (- 
0.0419 (- 
0.0396 (- 
0.0378 (- 



3-2^2 

-0%) 

-0.1%) 

-0.9%) 

-8.3%) 

15.4%) 

21.6%) 

26.7%) 

30.8%) 

33.9%) 



2. Null geodesies 

We now consider null circular geodesies, labeled by a 
"c" subscript. It is easy to show that, setting (5i = in 
equation (|5.4p . the dilaton field, H{r), gives no contribu- 
tion. Requiring Vc — V^ = Q we find 



2fc = rj:,. (5.9) 



The condition for a circular orbit reads 

V:'^2 ^'^f '"\ 2f,~rJ'J) , (5.10) 

whose sign does not depend on H{r). The orbital angular 
velocity is 




o _ 1 _ _ 

be t Tc 



(5.11) 



where he = L/E is the impact parameter. In the 
Schwarzschild case the condition ()5.9p restricts the null 
circular orbits to rc = 3M. The orbital angular velocity 
is MVLc = l/(3\/3). Results for a EDGB spherically sym- 
metric BH are summarized in Table HI In this case there 
is no effect from the coupling h and the difference between 
EDGB and GR theory is of the order 1%. From (jgl^ we 
find that V^' > and therefore only unstable null circular 
orbits exist. Hence also in EDGB case null geodesies are 
always unstable aganist small perturbations. 



TABLE II: Real and imaginary part of the QNM frequen- 
cies defined by the formula ujqnm — i^c I ~ i [n + 1/2)\\\ in 
the eikonal approximation [3J]. Terms between parenthesis 
are the differences with respect to the Schwarzschild case, for 
which MQ = MX = 1/(3^3). 





D/M 


MQc (±%) 


MX (±%) 


3 X 10"'' 


10"'' 


0.1925 (~ 0%) 


0.1925 (~ 0%) 


3 X 10"^ 


2 X 10"^ 


0.1925 (~ 0%) 


0.1925 (~ 0%) 


0.027 


0.013 


0.1925 (~ 0%) 


0.1925 (~ 0%) 


0.238 


0.129 


0.1926 (+0.1%) 


0.1924 (-0.1%) 


0.417 


0.246 


0.1930 (+0.3%) 


0.1921 (-0.2%) 


0.545 


0.350 


0.1936 (+0.6%) 


0.1916 (-0.4%) 


0.630 


0.441 


0.1942 (+0.9%) 


0.1909 (-0.8%) 


0.677 


0.516 


0.1948 (+1.2%) 


0.1901 (-1.3%) 


0.691 


0.572 


0.1953 (+1.5%) 


0.1892 (-1.7%) 



QNMs in the large I limit for spherica 
EDGB BHs 



symmetric 



Particularly promising for gravitational-wave detec- 
tion are the characteristic vibration modes of black 
holes [21I, [11]. These modes, called quasinormal modes 
(QNMs), are exponentially damped sinusoids and carry 
an imprint of the black hole, its features being indepen- 
dent on what exactly excited the modes. QNMs are ex- 
cited to a large amplitude for instance in the inspiral and 
subsequent merger of black hole or neutron star binaries. 
In fact, the ringdown phase of supermassive black holes 
can have a larger signal-to-noise ratio than any other 
signal m, [53, HI], justifying dedicated searches for 
ringdown in current gravitational-wave detectors such as 
LIGO or TAMA 59, 60]. 

It was recently shown ^] that there is a simple rela- 
tion between QNMs and circular null geodesies in general 
spacetimes in the eikonal, i.e. large-Z limit. 



iOQNM ^^lcl-i{n+l/2)\X\, (5.12) 
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where flc is defined in eq. (|5.1ip and 



A 



1.(2/. 



(5.13) 



is the Lyapunov exponent, describing the inverse of the 
instabihty timescale of the geodesic. For circular null 
geodesies the argument of the square root in eq. (|5.13[) 
is always positive. QNMs are easly computed from the 
equation above and again the difference between the real 
and imaginary part of the QNM frequency for a DBH and 
for a Schwarzschild BH is order 1%, as shown in Table 
im We have not attempted to compute the least damped 
QNMs for these black holes, but we will assume they 
behave in the same way as the coupling constant varies 
[3^. Discriminating such small percentage differences is 
in principle doable with the gravitational-wave detector 
LISA [21, 58], though it may be very challenging due to 
systematic errors [2l|, [H, [H, 113, HI]. As far as we know, 
the QNMs of DBHs have not been considered in the lit- 
erature, even though their computation is extremely rel- 
evant. Nevertheless, the computation of QNMs of black 
holes in higher-dimensional Gauss-Bonnet theories with- 
out a dilaton, or purely dilatonic black holes has been 
done [giI, [g^, [1^ , we expect the four-dimensional case to 
proceed along similar lines. 



B. Geodesies: rotating case 



TABLE III: Some particular cases of interest of the metric 
(I5.14P along the equatorial plane: axially and spherically sym- 
metric spacetimes. Then we specialize to the Kerr metric, and 
its slow rotation limit (Slow Kerr, SK) discussed in Appendix 
iBl and finally to the general slowly rotating metric described 
in Section ITVl and then we compare this to the slowly rotating 
Kerr solution, discussed. 



Axial Sym 


Sph. Sym 


Kerr 


SK 


Slow rot 


fir) 


fir) 


1 _ 'ML 


r 


fir) 


air) 


air) 


1_ 2M + £^ 


\ _ Ml 

r 


air) 


h{r) 




r 


^2 


^2 


j{r) 





2aM 
r 


2J_ 

r 




H{r) 


H{r) 








H{r) 



where again 61 = 1, for time- like and null geodesies 
respectively. Using the integrals of motion E and L the 
equations of motion read 



i 



V{r) = 5(r)e-2«('') x 

'h{r) - fir) - 2i{r)EL 
j2(r)-f /(r)%) 
3ir)E + f{r)L _^(,) 
f{r)+fir)h{r) 
h{r)E-jir)L _^(,) 
fir) + fir)hir)'' 



(5.17) 



We will now consider geodesies in the the slowly rotat- 
ing DBHs spacetime studied earlier in Section IIVI For 
that, let's take the general four-dimensional static ax- 
isymmetric spacetime on the equatorial plane 

= fir)dt^ - ^dr^ - h{r)dip^ + 2jir)dtd(p , (5.14) 

gir) 

where we specialize to 9 — 7r/2 and d9 — 0. We also 
consider the dilaton field on the equatorial plane and set 
H{r) = 26(/)(r). Some interesting special cases of the 
metric above are listed in Table IIIII Following Chan- 
drasekhar (5^ we consider planar orbits, for which the 
Lagrangian is 

2C = e"('-^ (fir)P - _ j^^^^^2 ^ 2j{r)i^\ . 

V din J 

(5.15) 

The generalized momenta are 

e^W [f{r)i + 3{r)^]^E, 
e^W [-j{r)i + hir)^]^L, 



Pt 


— const — 


-Pv 


= const = 






-Pr 


— 

gir) 



and the Hamiltonian is 



2n = 2{pti + p^i^ + pri- - C) ^ Ei-L(p- 



e"^r) 

air) 



(5.16) 



Time-like geodesies 



Setting (5i = 1 in eq. (|5.17p and requiring V = V' = 
at the radius r — rt we find a system of non-linear 
algebraic equations for E and L 



ht E' 



e"*if^+ftht) 



ftL^-2jtEL 
f[L^~2j[EL + 
e"' [2 Jt ft + ft ht + ftK + H[ (jf + ft ht)] 



The system above can be solved analytically but the form 
of the solutions is not particularly useful. Solutions of the 
system above can be substituted in 



fi'L^~2f'EL~e 



Ht 



ft+ftht 
X {[H'^ + {H[mt+ftht)+ 
+2H[i2jtf + f'tht + ftK) + 

+2 jf -t- 2 Jt ft' + ft ht + 2 /; h'^ + ft h'^))] (5.18) 

and in 

if ^ jt E/L + 
i ht E/L - jt ' 



From eq. (|5.18p above we can find the ISCO for a generic 
rotating BH asking for V/' < 0. Setting /i(r) = and 
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TABLE IV: Results for co-rotating time-like (at the ISCO) and null geodesies in slowly rotating EDGB BH spacetimes compared 
to the slowly rotating Kerr case. 



q' /M^ 










3 X 10~^ 
3 X 10"^ 
3 X 10"^ 


0.3 
0.2 
0.1 


0.0901 (~ 0%) 
0.0808 (~ 0%) 
0.0737 (~ 0%) 


0.0900 (-0.1%) 
0.0807 (-0.1%) 
0.0736 (-0.1%) 


0.2220 (~ 0%) 
0.2101 (~ 0%) 
0.2005 (~ 0%) 


0.027 


0.3 


0.0901 (~ 0%) 


0.0896 (-0.6%) 


0.2220 (~ 0%) 


0.238 


0.3 


0.0903 (+0.3%) 


0.0855 (-5.1%) 


0.2224 (+0.2%) 


0.417 
0.417 


0.3 
0.1 


0.0910 (-fl.0%) 
0.0742 (-fO.6%) 


0.0818 (-9.2%) 
0.0668 (-9.3%) 


0.2232 (+0.5%) 
0.2012 (+0.4%) 


0.545 
0.545 


0.3 
0.1 


0.0919 (-f2.0%) 
0.0746 (-H.2%) 


0.0787 (-12.6%) 
0.0641 (-13.0%) 


0.2244 (+1.1%) 
0.2019 (+0.7%) 


0.630 

U.DoU 

0.630 


0.3 

A O 

0.1 


0.0929 (-h3.2%) 
0.0751 1+1.9%) 


0.0763 (-15.3%) 

A Af;Q1 f 1 CC Q07'\ 

U.UDoi \ — 10.0/0) 

0.0619 (-16.1%) 


0.2258 (+1.7%) 

noioA ( 1 1 QOy\ 
U.zizy y+l.o/o) 

0.2027 (+1.1%) 


0.677 
0.677 
0.677 


0.3 
0.2 
0.1 


0.0939 (-f4.3%) 
0.0834 (-^-3.2%) 
0.0756 (+2.5%) 


0.0744 (-17.4%) 
0.0661 (-18.1%) 
0.0600 (-18.6%) 


0.2272 (+2.4%) 
0.2139 (+1.8%) 
0.2035 (+1.5%) 


0.691 
0.691 
0.691 


0.3 
0.2 
0.1 


0.0947 (+5.1%) 
0.0839 (+3.8%) 
0.0759 (+3.0%) 


0.0729 (-19.0%) 
0.0647 (-19.9%) 
0.0586 (-20.4%) 


0.2284 (+2.9%) 
0.2148 (+2.2%) 
0.2041 (+1.8%) 



j (r) = equations above reduce to ones for a spherically 
symmetric spacetime. The requirement for the energy 
E to be real imposes a constraint rj > tmin- Results 
for circular time-like geodesies are shown in Table [TVl for 



co-rotating orbits. For tmin < rt < nsco only unstable 
circular orbits are permitted. For rt > nsco circular 
orbits are stable, while for rt < ^min no circular orbits 
exist. 



Measurements of the ISCO can be used to evaluate the 
angular momentum of an astrophysical BH [28]. Figure 
[S] shows how the ISCO and the orbital frequency, fit de- 
pend on J for slowly rotating BHs both in EDGB theory 
and in GR. Figure [8] can be used to evaluate the angular 
momentum once the ISCO has been measured (see Fig. 
2 in Ref. j^Si] for details). Again the role of the coupling 
b is leading for timelike geodesies. As a qualitative re- 
sult the orbital frequency is smaller for a slowly rotating 
EDGB BH than the one for a Kerr BH and the ISCO is 
larger. Differences range from 10% to ^ 20% depending 
on the angular momentum J /M'^ and on the dilatonic 
charge D/M . The difference is monotonically decreasing 
for larger rotations, it strongly depends on h but not on 
J. The effect of rotation seems to be subleading. 



2. Null geodesies 



Focusing on circular orbits {5i — 0), we require V = 
y = at the radius r — rr- These two conditions read 



E 
T 




'.;(|)'-/;+2i;f.o^ 



In this case 



fc+fchc 



h'! 



E 



(5.19) 
(5.20) 

(5.21) 



which is positive at r = rc. As expected the positiveness 
of equation above does not depend on the dilaton H(r). 
Therefore only unstable null circular orbits are allowed. 
The angular velocity is 



J_ ^ _^ y/jg + fchc 
be hr. hr 



(5.22) 



where he = L/E is the impact parameter. The double 
sign in the above equation is related to orbits which are 
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FIG. 8: Top Panel: The innermost-stable-circular orbit ra- 
dius, risco/M as a function of the angular momentum, 
J/M^ for a slowly rotating DBH with q'/M^ ~ 0.691 for 
b — 1/2 (solid line), fe = (dashed line) and in the GR limit, 
a' /M^ ~ (dotted line). From these plots and from the 
measurements of risco the value of the angular momentum 
J/M^ can be evaluated (see Fig. 2 in Ref. [23|). Bottom 
Panel: the orbital frequency at the ISCO f2isco as a function 
of the angular momentum, J/M^ for a slowly rotating DBH 
with a'/M^ ~ 0.691 for b = 1/2 (solid line), 6 = (dashed 
line) and in the GR limit, a' /M^ ~ (dotted line). Negative 
values for J correspond to counter-rotating orbits. 

co-rotating and counter-rotating with the BH. Results 
for co-rotating null orbits are shown in Table IIVI and in 
Fig. [HI As we expect, results for J/M^ smoothly 
tend to the ones for the non-rotating case. Interestingly 
enough, percentual differences between GR results can be 
larger than order 5% for co- rotating orbits. The radius 
for circular null orbits is always smaller for a rotating 
EDGE BH than for a Kerr BH in GR. 



VI. DISCUSSION 

Einstein-Dilatonic-Gauss-Bonnet theories are viable 
theories of gravity, which share many features in com- 
mon with Einstein's gravity but have a better under- 
stood quantum limit. We have improved the linear sta- 




-0.3 -0.2 -0.1 0.0 0.1 0.2 0.3 



FIG. 9: The orbital frequency, Qc, for circular null-like orbits 
as a function of the angular momentum, J/M^ for a slowly 
rotating DBH with a'/M'^ ~ 0.691 (solid line) and in the 
GR limit, a'/A/^ ~ 0.03 (dashed line). Negative values for J 
correspond to counter-rotating orbits. 

bility analysis for static black holes in this theory. We 
found that they are stable and can in principle be used 
to discriminate between the two theories. For the slowly 
rotating black holes we studied here, the differences in 
measurable quantities amount to a few percent in the 
most favorable cases (large coupling constant). Given 
the current state-of-the-art, it does not seem possible to 
distinguish the correct theory of gravity from measure- 
ments of either the ISCO (from EM observation) or even 
of ringdown modes with future gravitational wave detec- 
tors. 

There are at least two important extensions to be 
made, which could prove to be very relevant: 

(i) Consider highly spinning black holes. One could 
attempt to extend the small rotation expansion to second 
order in Q, but we found such procedure to be extremely 
complex, due to the symbolic manipulations required. If 
carried through, it would allow a computation of other 
multipole moments of the hole, which could potentially 
lead to important tests [2l|, Given that many of 
the astrophysical black hole candidate have a large spin, 
this could be a promising way of discriminating EDGB 
and Einstein's theory. To find a highly spinning DBH 
solution most likely requires a full numerical integration 
of the field equations (see for instance Ref. [66| for a 
concise overview of possible methods to handle this). 

(ii) Waveforms of inspiralling particles. The calcula- 
tion of waveforms is highly non-trivial: it includes either 
a re-derivation of the PN expansion, or a full-blown nu- 
merical solution of the problem. The pay-off is huge: 
particles orbiting around black holes are able to probe 
the background geometry in its entirety (67j . The result- 
ing gravitational- waveform should be sensitive enough to 
the field equations, specially at late stages in the inspi- 
ral, when the particle is about to merge. Indeed, the 
particle spends a large number of cycles near the ISCO, 
which could potentially increase by orders of magnitude 
the ability to probe the geometry better, if compared to 
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ringdown. 

Of course much more remains to be done. For instance, 
it would be desirable to have a quantitative analysis of 
the lowest lying ringdown modes. Most of all, it would 
be extremely important to understand how general the- 
ories of gravity affect the strong field regime, the ISCO 
location and frequency, etc. Perhaps then one might un- 
derstand how to discriminate between the correct theory 
of gravity with future electromagnetic or gravitational- 
wave observations. We hope to tackle these issues in 
future works. 
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APPENDIX A: LINEAR STABILITY ANALYSIS: 
AXIAL PERTURBATIONS 

In this appendix we derive the Schrodinger-like equa- 
tion which represents a eigenvalue problem for the com- 
plex frequency lo. Detailed calculations can be found in 
the online MATHEMATICA notebooks which are avail- 
able at http : //paolo . casadiale . com/EDGB_BHs . zip. 

The coefficients in eqns. p.3p - p.5p read 



i e-^+^ (2e^ (A' - V) - a'e'^° {V {-cj^p' {-3A' + 2cl,o' + T) - 2(/.o") - 2<j>o'r")) 



uj{2e^ + a'e't'a 



A 



'4>o - 2 {4>a'^ + 0o" 



2ruj (e^r - a' e'>'» (j)o') 



2e^r(l)'n^T' - 2e^ rP 



+ (2e^ + 1) - 2) + r (r^J,' + (2 + rF') (F' - A') + 2rF"))] } . 



-2eW ' -r^'^"^ 



r 



(Al) 
(A2) 



Vo" + 0o' (2e^rLu^ - ((e^ (1(1 + 1) - 2) - 3rA') F' + rP'' + 2rF")) 



(A3) 



The explicit for of the Cij coefficients is not shown here. 
One can numerically prove that equation (13. 5p is au- 
tomatically satisfied as a consequence of the other two 
equations and of background solutions. So we are left 
with a system of two ODEs for the unknown functions 
hpi^r) and hi{r). Eliminating hp from the first order equa- 
tion (|3.3p we obtain a second order differential equation 
for hi, 

A{r)h'l{r) + 2Bir)h[{r)+Cir)hi{r) =0, (A4) 

where A{r) ^ A2, B{r) = + A'^ - ^) and C(r) = 
~ ^7^ ~ -^1,;- The function C(r) can be decomposed 
in C(r) = Q{r) +u;'^E{r), where Q{r) and E{r) don't 
depend on uj. The functions A,B,E and Q depend on 
the radial background functions F, F', F", A, A', (j>o, 0q, (pQ 
and on ^q". 

In order to eliminate the term proportional to h'l we de- 
fine 



F — exp 



(A5) 



Setting u = Fhi, equation (|A4p takes the Schrodinger- 
like form 

u"{r) + [V{r) + uj^K{r)] u{r) = , (A6) 
which is eq. (|3.6p with 



E 
A 



In the limit r the coefficients A, B, Q and E take 

the form 









(1 







A = A,, + 0{r-rh), B 



By, 



Q 

E 



{r - rh)^ 
Eh 1 



O 



1 



r - rii 
O 



1 



0(1),(A8) 
(A9) 
(AlO) 



71 {r-Th)^ 

In equations above we have used the following asymptotic 
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behavior near the event horizon 



"AM ^ 
0(r) = 



Ai(r - Th) + ... , 
li{r - Th) + ... , 
0ft + <P'h{f - rh) 



(All) 
(A12) 
(A13) 



where F is the hypergeometric function and G is the 
Meijer function. The asymptotic behavior of equation 
(|B1|) is the same as eq. (j4.6|) . The asymptotic behavior 
(|4.5p and the solution at infinity (|4.7p imply that only 
the I — 1 mode is allowed. A list of the general solution 
of the equation above is: 



where 



Ai - 



(A14) 



and 7i is an arbitrary finite positive integration constant, 
which cannot be fixed by the equations of motion, since 
the latter involve only r'(r) and not r(r). This constant 
is fixed by the asymptotic limit of the solutions at infinity. 
The asymptotic behavior near the event horizon for V{r) 
and K{r) is 



V 



Vh 



K 



(r - ThY 



where Vh ^ 1/4 and Kh < -iA'P. The constants 
Ah, Bh,Qh, Eh, Vh and Kh depend on r^, (j)h and ^J^. For 
a Schwarzschild BH, Vi = 1/4 and Kh = ^AP. When 
r — !■ oo, — > and K ^ 1. The asymptotic behavior 
for equation p.6p is described by Eqs. p.7p and (|3.8p . 



APPENDIX B: SLOWLY ROTATING KERR 
BLACK HOLES IN THE HARTLE 
APPROXIMATION 

A useful check on the method described in Section IIVI 
is to consider classical GR, and repeat the calculation 
adding rotation to a Schwarzschild black hole. Perturb- 
ing the Schwarzschild metric and retaining only first or- 
der terms in Q, equation (|4.3p reduces to 

n['{r) + ^n[ + ^-^^^ + % , = o, (Bi) 

r r(r — 2M) 

which can be solved analytically in terms of special func- 
tions 



n,M = Ci F(l-i,2+i,4;— j+Cj GiS 



(- 



-1-1,1 
-3,0 
(B2) 



^i{r)=Ci + ^, 1 = 1. 



Co 



2M 



I = ; 



Af2 



3M 



^^3M=C,(^-2)(3^-4 
M 



(B3) 
3 

2 + 
(B4) 



AP 

4^ + 10^ 



45 / r 4 
T \M ~ 3 



30 105-1-45 — 

r AI 



\M J 



log 



r - 2AI 



,1 = 3(B5) 



Demanding a correct asymptotic behavior at infinity, 
eq. (|4.5p . we verify that / = 1 is the only mode able 
to satisfy regularity at the horizon (the result extends to 
I > 3 modes, though we do not show those here). We 
therefore get the exact result = 2J/r^. This coincides 
with the expansion (first order in a = J/M) of a Kerr 
metric in Boyer-Lindquist coordinates. 



ds' = 



2M\ , o / 2Af\ ^ , ^ 
d9^ + sin^ e (dip - ^dtf 



(B6) 



which we have named Slow Kerr (SK) metric in this work. 
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